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Abstract 

The radiative spin-flip transition rates of heavy quarkonium states 
depend sensitively on the matrix elements of the effective confining 
interaction through the associated two-quark exchange current opera¬ 
tor. The Hamiltonian model based on a scalar linear confining interac¬ 
tion with a single gluon exchange hyperfine term is shown to provide 
an adequate description of the J/fj —> r]cl and —> yd decay 

widths once the relativistic single quark magnetic moment operator is 
treated without approximation and the exchange current is taken into 
account. Predictions are given for the radiative spin-flip decay widths 
of the IS*,25 and 35 states of the cc, bh and Bf systems. In the Bf 
system the gluon exchange current also contributes to the spin-flip 
transition rates. 
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1 Introduction 


The spectra of the cc and bb meson systems are fairly well reproduced by 
a Schrodinger equation based description with an interaction Hamiltonian 
formed of a scalar linear conhning interaction and a single gluon exchange 
model for the hyperhne interaction il- This model is qualitatively sup¬ 
ported by numerical construction of the effective interaction by lattice meth¬ 
ods Bi- For this linear conhnement -|- gluon exchange model the parameter 
freedom has recently been narrowed considerably by numerical precision de¬ 
termination of the effective quark-gluon coupling strength p, |^. 

While the spectra of the cc and bb systems are fairly well accounted for by 
the simple linear conhnement -|- gluon exchange interaction model, the situ¬ 
ation concerning the radiative widths of the spin-hip transitions J/V' —> rjc'y 
and 'ip{2S) has remained unsettled i I, I- We show here that these 

radiative transitions strengths are in fact also satisfactorily described by the 
conventional quantum mechanical framework, provided that the full Dirac 
structure of the quark current operators is taken into account along with 
the exchange current operator that is generated by coupling to intermediate 
negative energy states by the linear scalar conhning operator Fig.|^, once the 
same Hamiltonian is used to generate both the wave functions and the cur¬ 
rent operators. This exchange current contribution is decisive for achieving 
agreement between the calculated and measured Ml transition rates of the 
J/'ip and the ip' . The exchange current operator is required by current con¬ 
servation with the conhning interaction. 


The contribution of the exchange current operator that is associated with 
the linear potential is essential for reproduction of the empirical decay widths. 
Moreover this operator only arises under the assumption that the ehective 
conhning interaction is a Lorentz scalar (“S'”), which supports the consis¬ 
tency requirement of ref. |^. If the Lorentz invariant structure of the con¬ 
hning interaction were a vector invariant (“H”), the associated exchange 
current operator would contain no spin-hip term when the quark and anti¬ 
quarks have equal mass as in the cc system, and agreement with the empirical 
decay widths would be excluded. 
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Figure 1: Exchange current operators associated with the effective scalar 
confining and gluon exchange interactions, with intermediate virtual qq ex¬ 
citations. 


In the case of the B 
0 


observed 


system, the first state of which has recently been 
the exchange current that is associated with the vector gluon 


exchange interaction in contrast also contains a component that is antisym¬ 
metric in the spins of the quarks and antiquarks, which in this case have 
unequal mass. That exchange current consequently also plays a role in the 
spin-flip transitions of the system, because of the large difference in mass 
between the c and b quarks. The role of these exchange currents is shown 
here by a calculation of the the decay widths for the radiative spin-flip tran¬ 
sitions B^{J = 1, nS) —>• B^{J = 0, 1S')7 {n = 1,2,3). The decay widths 
are calculated here with full account of both the confinement and gluon 
exchange currents. The radiative decay widths have been calculated in the 
non-relativistic impulse approximation without consideration of the exchange 
current contributions in ref.|]T3. Experimental determination of these decay 


widths should be decisive for settling the quantitative importance of the ex¬ 
change current operators and the validity of the present phenomenological 
approach to the structure of heavy quarkonia. 
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The B'^ systems differ from the cc and bb systems also in that the 
meson have magnetic moments. These also obtain considerable exchange 
cnrrent contribntions which are calculated here. 


The framework used here is similar to that employed in refs. Bi for the 
calculation of the rates of the Ml transitions in the cc and bb systems. The 
affirmative conclusion about the signihcance of the exchange current con¬ 
tributions here is new. The numerical values differ from those in refs.[^ 
partly because the unapproximated Dirac magnetic moment operator is con¬ 
sidered here, and because of an overestimate of the lowest order relativistic 
correction to the quark magnetic moment operator. 


This paper is divided into 5 sections. In section § we describe the Hamil¬ 
tonian model for the heavy quarkonium systems, and demonstrate that it 
describes the spectra of the cc and bb mesons satisfactorily. In section |0 
we show the predicted spectrum of the B^ system. In section ^ we describe 
the model for the current operator for the heavy quarkonium system. The 
corresponding calculated decay widths for the radiative spin-flip transitions 
of the cc, bb and B'^ mesons are presented in section Finally section | 
contains a concluding discussion. 


2 Model Hamiltonian 

2.1 Hamiltonian 

The linear scalar conhnement -I- gluon exchange hyperhne interaction Hamil¬ 
tonian model for heavy quarkonia Qq, is formed of quarks Q and antiquarks 
q with in general unequal mass 

H = Hkin + Vc + VG ( 1 ) 

where H^in is the kinetic energy term and Vc and Vq are the potentials that 
describe the conhning and gluon exchange interactions respectively. 

The kinetic energy term, which in heavy quark effective theory is taken 
to include the terms to order is 
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( 2 ) 


Hkin = mq + rriq + 


p 


2mr 


1 

8 


V ^ ’ 


where is the reduced mass = mQrriq /{uiq + rriq) and p is the relative 
momentum. The term of Cl(p^), which will be treated as a perturbation in¬ 
creases in importance with excitation number, and is non-negligible because 
of the compact spatial extent of the conhned quarkonium wave functions. 


The scalar conhning potential has the form (to order m ^): 
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cr 
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The last antisymmetric spin-orbit term in this interaction vanishes for equal 
mass quarkonia. 


The single gluon exchange interaction potential to order m ^ has the form 
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mqmq 

Here Su = 3aq ■ faq ■ f — aq ■ is the usual tensor interaction operator. 


The interaction operator constructed to order by lattice methods in 
refs. P, ^ contains in addition to the operators of the form contained in the 
model Hamiltonian O that involve the non-local operator p/m to second 
order, but which do not signihcantly affect the calculated spectra. We do 
not include such here, mainly because they lack an immediate dynamical 
interpretation. The interaction operator derived in refs, p, ^ has a string 
tension c of ~ 1.1 GeV/fm, with a slight flavor dependence. 
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cc 


bb 

c 

1.11 GeV/fm 

1.11 GeV/fm 

1.11 GeV/fm 

(Ts 

0.40 

0.32 

0.29 

me 

1.48 GeV 

1.48 GeV 

- 

mb 

- 

4.72 GeV 

4.72 GeV 


Table 1: Parameter values for the model Hamiltonian (|^ used here. 

We shall here use parameter values close to those suggested in 
but determine the precise values by phenomenological hts to the cc and bb 
spectra. The values used are listed in Table |^. The values for the color 
hue structure constant ag in the table are close to those obtained by lattice 
methods in refs. [^, |^: ag = 0.38 at the charmonium scale rric ~ 1.3 and 
ag = 0.22 at the bottom scale mf, ~ 4.1 GeV respectively. These values agree 
well with the inverse logarithmic momentum dependence predicted by per¬ 
turbative QCD. The phenomenologically determined quark masses in Table 
in should be viewed as constituent quark masses. The overall description of 
the empirical spectra corresponds to that achieved in refs, p, Q. 


2.2 The cc and bb Spectra 


The calculated spectra for the cc and bb systems that are obtained by solving 
the Schrodinger equation for the model Hamiltonian (m) are shown in Figs. 
^ and ^ respectively. In the calculation the term of order in the kinetic 
energy operator (Q), as well as all the terms of order in the interaction 
potentials W (0) and Vq (|) were treated in hrst order perturbation theory. 
The parameter values used were those in Table |l|, which are close to those 
employed in ref. [^. The perturbative treatment is motivated by the fact 
that the main focus here will be the calculation of the radiative widths of the 
F-states, for which the only terms in the hne structure part of the interaction 
that matter are the delta function terms in and (^), which have to be 
treated perturbatively. 
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state 

s = 0 

J = l 

j = 1-1 


j = 1 + 1 

IS 

2978 


3159 



(2979) 


(3097) 


2S 

3558 


3683 



(3594) 


(3686) 


3S 

3895 


4004 


4S 

4130 


4231 


IP 

3547 

3471 

3538 

3567 


(3526) 

(3415) 

(3511) 

(3556) 

2P 

3903 

3827 

3892 

3926 

ID 

3796 

3798 

3802 

3790 

2D 

4079 

4075 

4083 

4077 


Table 2: The cc-states in MeV. 


Figure 2: Charmonium spectrum with nic = 1480 MeV obtained from the 
Hamiltonian Charmonium Spectrum (State energies in MeV) 
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state 

s = 0 

3 = 1 

j = 1-1 

s = 1 

3 = 1 

j = 1 + 1 

IS 

9411 


9490 





(9460) 


2S 

9964 


10008 





(10023) 


3S 

10311 


10346 





(10345) 


4S 

10589 


10621 


IP 

9895 

9868 

9890 

9903 



(9860) 

(9892) 

(9913) 

2P 

10248 

10224 

10243 

10256 



(10232) 

(10255) 

(10268) 

ID 

10138 

10135 

10138 

10140 

2D 

10435 

10431 

10435 

10438 


Table 3: The 66-states in MeV. 


The quality of the agreement between the calculated and empirical cc and 
66 spectra in Figs. ^ and |] is similar to that of the corresponding spectra ob¬ 
tained by the full numerically constructed Hamiltonian model [^, In both 
cases the splitting between the lowest S'-states {rjc, J/ip) is overpredicted and 
that of the P-states {Xcj) is somewhat underpredicted. A better description 
of the rjc — J/'ll) splitting may be obtained with the Buchmiiller-Tye potential 
used in ref. P, but only at the price of underprediction of the excited 


states. This problem appears to be generic to Hamiltonian models of the type 
O- The calculated spectra are nevertheless in good overall agreement with 
the empirical spectra, with hyperhne splittings of the multiplets which are 
ordered as the empirical ones. The model therefore appears realistic enough 
for a quantitative calculation of the spectrum of the P+ system as well . 


In Fig. 1^ the empirical cc state at 4160 MeV has been given the assign¬ 
ment 4^S'! state, although it may be a 2D state [113. Similarly the empirical 


T state at 10580 MeV has been given the assignment 4^S'! in Fig. ^ These 
states fall sufficiently low for this assignment with the present model, mainly 
as a consequence of the term in the effective conhning interaction (^). 




















The Bottomonium Spectrum (State energies in MeV) 
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Figure 3: Bottomonium spectrum with rrih = 4720 MeV obtained from the 
Hamiltonian (1) 
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2.3 The Meson Spectrum 

The first state of the B'^ meson spectrum, formed of c and b quarks - pre¬ 


sumably a I'^S'i-state - has recently been observed [111. The Hamiltonian 


model (]I|) may be applied directly to the calculation of the spectrum of the 
system. By setting = 0.32 the IS'-states agree with the earlier spectra 
calculated in ref. P|. The precise value for the empirical mass of the discov¬ 
ered H^-state is not yet known. 


The predicted meson spectrum is shown in Fig. This overall fea¬ 
ture of this spectrum is that of an interpolation between the spectra of the 
cc and the bb systems. The large mass difference between the c and b quarks 
brings no qualitatively new features, although it does lead to the antisym¬ 
metric spin-orbit interactions in (j^) and (^. 
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The predicted energy values in Table ^ and in Fig. agree fairly well with 
those obtained in ref. |^, although the major shell spacings here are some¬ 
what smaller, a feature that was built in because that is in better agreement 
with the empirical cc and bb spectra. 


3 The Current Operator 

3.1 Single Quark Current 

Under the assumption that the b and c quarks are point Dirac particles, their 
current density operators have the form 


< p'\j{0)\p>=iQeu{p')^u{p), (5) 

where Q is the quark charge factor (-1-2/3, —1/3 respectively). Assuming 
canonical boosts for the spinors u, it, the spin part of the current yields 


the following magnetic moment operator |T5 


p = Q 


rrir, 


mq) A/l + n 




1 -^ 1 - 


y/TTW, 


Pn- 


( 6 ) 


Here v is the quark velocity v = {p' -|- p)/2mQ, rup is the proton mass and 
Pat is the nuclear magneton. 

For the cc and bb systems the spin-magnetic moments (§) combine to a 
spin-flip term of the form 


= 


2 1 
3’~3 


rrir 


_ ~ 

mq) 


1 - 


1 - 


vTT^. 


Pn- 


(7) 


The matrix element of this operator for a spin-flip transition J = 1 J = 0 
between F-states of the cc and 66 system may be expressed as 


< /l/7i|i >= 4 


roo , roo 

dr'r ^ / 

3 Jo 

q 


3’ 3 

drr‘^Pf{r') 


rrip 

1 


Pn 


dpp^ 


dz 


'-1 


a/1 + 






vT-jT" 


•r?2 


Jo r 


p+ 


Jo r 


p 




( 8 ) 
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Here v = p/m, and (pi{r) and ^Pf(r') are the radial wave functions for the 
initial and hnal S'-states. 

The matrix element (§) would be considerably simplihed by expansion 
to lowest order in v of the square roots in the expression (D. This would 
however lead to misleading results as the radiative widths of the S'-states are 
very sensitive to the model of the current operator and the wave functions 
1^). For reference we note that if only the order (r;)° term is kept in (|^ the 
expression (||) reduces to the standard non-relativistic impulse approximation 
result 


(|-i) *"V/(rWo(f)v>.(r). (9) 

Here Q = c, b for cc and bb respectively. The higher order terms in (|^) lead 
to numerical values that differ signihcantly from this result as shown below. 


If the magnetic moment expression is expanded to second order in 
the correction factor to the non-relativistic magnetic moment operator is 
1 — 2r;^/3. The coefficient 2/3 in this bracket was obtained incorrectly as 5/6 
in refs. Bi. with a consequent overestimate of this relativistic correction. 
The factor 2/3 is obtained as the sum of two terms, one of which arises from 
the normalization of the Dirac spinors (1/2) [|T^ and another, which arises 
from the difference between the initial and final quark momenta (1/6). 


For the system the spin dependent part of the magnetic moment 
operator (^) is modified to 


MBc) = Qmp{ac + ai)< 


1 -- 1 






1 


1 -- 1 - 


>Pn 


\/l + T 


c, 


+ ^mp{ac - ai)- 


^c^l + v^ 


T12 

C 


1 -^ 1 - 




.-r2 
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1-^1 


+ T 


b . 


\/l + T 

f^N, 


Tt2 
c , 


( 10 ) 
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where Vc = p/mc and hfe = p/mh respectively. The term that is symmetric 
in the spins gives rise to a magnetic moment for the ^Si state of the 
system, whereas the term that is antisymmetric in the spins gives rise to 
spin-flip transitions. 

The radiative widths calculated using these single quark magnetic mo¬ 
ments - the impulse approximation - do not agree with the presently known 
empirical rates for Ml transitions. These discrepancies can be considerably 
ameliorated, if not entirely eliminated, by taking into account the exchange 
current contributions. 

3.2 Exchange Current Operators 

The scalar conhning interaction may excite virtual quark-antiquark states, 
which are deexcited by the external electromagnetic held. This process, 
which is illustrated by the Feynman diagrams in Fig. |I], generates an ex¬ 
change current - or two-quark pair excitation current, which in the cc and 
bb system may, to order be expressed in compact form as 

CT 

UC) = -ji, (11) 

mg 

where ji is the corresponding single quark current operator. This exchange 
current operator is required to satisfy the continuity equation once the terms 
of order 1 /m^ are included in the conhning interaction (^. 

The spin dependent part of the corresponding exchange magnetic moment 
operator is then, to lowest order in 1 /m^, for the cc and bh systems: 

P2{C) = - - {aQ-aq)pN, 12 

V3 3/ tuq mq 

where the factor 2/3 in the bracket applies to the cc system and the factor 
— 1/3 to the bb system and mq = nic and mh respectively. The matrix element 
of this operator for a transition z —^ / -|- 7 is then 

< f\MC)\i >= - drr^pf{r)crjo ipi{r). (13) 

The contributions of these matrix elements turn out to be essential for the 
explanation of the empirical decay widths of the transitions J/z/ —>■ rjpy and 
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In the case of the system the proportionality between the conhnement 
exchange current and the single quark current operator (^) is lost. The 
complete scalar pair excitation current operator for the B'^ system, to order 
has the expression 

erf pc + pc' Pi + Pi 

The spin dependent part of the corresponding magnetic moment operator 
takes the form 


cr 


fl2{C)[Bt] = -IT 


rrir. 


2 1 
—H-^ 


rat 


rat 


(<?c + ^h) + 


, rat 


rat 


(<?c — <?&) \ ^N- 


(15) 


The hrst term in this expression contributes to the magnetic moment of the 
J > 0 B'^ mesons. 

In the B'^ system the gluon exchange current operator that arises from 
excitations of virtual qq pair states will contribute to both the magnetic 
moments and to the spin-flip transition strengths. The complete expression 
(to order of this gluon exchange current operator is (in momentum 

space): 


UG)[Bt] = e 

^ 1 Pc + Pc' 

3 /c^ mcTah 


87ra, f 2 1 pi+pi 
9 \ 3 rnemb 

+ i + 

\mcmb ml J 


+ i 


X kr 


CTr 


+ 


(Tb 


X kh 


mt mcmb J 


(16) 


Here the momentum operators ki and kc denote the fractional momenta de¬ 
livered to the h and c quarks respectively {q = ki + k^). 


The spin dependent part of the magnetic moment operator for the gluon 
exchange current operator ([I^) takes the form 


-2 I + ^h) 


MG)[Bt] = ^M {—,+ 

27 r ( \mt. 


mcmb ml 
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( 17 ) 


Note that the presence of the spin-flip term in (|^ is solely a consequence 
of the difference in mass of the charm and beauty quarks. 


4 Radiative Spin-Flip Transitions 


The spin-flip part of the magnetic moment operators considered above may 
be expressed in the general form 


hS.F — “ ^q), 


(18) 


where Ai is the matrix element of the orbital part of the operator for the 
transition and e is the elementary charge. The width for a radiative spin-flip 
transition of the form Qq{J = 1) —>■ Qq{J = 0)7 then takes the form 



(19) 


Here ctem is the hue structure constant, q is the photon momentum in the 
laboratory frame, and Mi and Mj are the masses of the initial and hnal states 
respectively. The matrix element M. is formed of a single quark current term 
and an exchange current term. These expressions for the heavy quarkonium 
systems are given explicitly below. 


4.1 Spin-flip transitions in the cc system 

For radiative spin-flip decays of the form '0(nS') —>• 7c7 and r]'^{nS) -07 the 

matrix element (|^ M. may be derived from the single quark operator ( 0 ) 
and the exchange current contribution ([I^) that is associated with the scalar 
conhning interaction. The matrix element M. then takes the form [0: 



( 20 ) 


The factor l/SrUc arises from the charge factor 2/3 and the factor l/2mc in 
the magnetic moment operator. The dimensionless integrals Ji and Ic that 
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arise from the single quark operator (|^) and the conhning exchange current 
([T^) respectively have the explicit expressions (cf.(|[), (p!3[) ). 

POO p1 poo poo 

Ji = 4 / dpp^ / dz dr'r'’^ / drr‘^ip*Ar') 

Jo J-i Jo Jo 

VTTt^ { 3 ( a/1 + 

jo (^r'^Jp‘^ + pqz/2 + q^/16^ jo (^rp‘^ - pqz/2 + q^/16^ pi{r), (21) 

Ic = -/ drr^(p*f{r)r(pi{r). (22) 

rric Jo 

Here v = p/m. The hrst of these two integrals reduces to the form (H) (with 
exception of the factor 2mpp/3mQ) in the static limit n —>■ 0. 



Transition 

NRIA 

RIA 

RIA-f-conf 

Exp. 


2.87 keV 

2.18 keV 

1.05 keV 

1.14 ± 0.39 keV 

2^5i ^ D5o 
23^1 ^ 21^0 

1.54 keV 
1.43 keV 

0.206 keV 
1.03 keV 

1.32 keV 
0.14 keV 

0.78 ± 0.24 keV 

2^50 ^ 

0.188 keV 

0.469 keV 

0.351 keV 

— 

3351 ^ 21^0 
33^1 ^ 3^So 

0.795 keV 
0.61 keV 
2.32 keV 

0.172 keV 
0.0753 keV 
1.57 keV 

0.929 keV 
1.24 keV 
0.0124 keV 

— 

3^50 ^ 

31 Ao ^ 23^1 

0.167 keV 
0.0133 keV 

0.228 keV 
0.12 keV 

0.372 keV 
0.233 keV 

_ 


Table 5: Decay widths for the cc system. The columns NRIA and RIA con¬ 
tain the results of the non-relativistic and relativistic impulse approximations 
respectively. The net calculated decay width obtained by combining the rel¬ 
ativistic impulse approximation and exchange current contributions is given 


in the column RIA-t-conf. The empirical values are from ref. 


The calculated decay widths for spin-flip transitions between the 5—states 
of the cc system are given in Table |^. In Table ^ the decay width results that 
are obtained in the non-relativistic and relativistic impulse approximations 
are also given both in order to facilitate comparison to the results of refs. 
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1^, 1^ and to bring ont the significant role of the exchange cnrrent contri- 
bntion. Note that in the nnmerical calcnlation the empirical valnes for the 
photon momenta have been nsed in order to obtain the correct phase space 
factors. This makes a difference only for transitions that involve the J/ijj 
state, which is somewhat overpredicted by the present model. This issne 
does not affect the wave fnnctions, which were calculated without account 
of the hyperhne terms of order which were treated in lowest order 

perturbation theory. 

The results in Table ^ show that the non-relativistic impulse approxi¬ 
mation overestimates the rate for J/-^ —> pc7 by a factor ~ 3 as has been 
noted earlier BB- By employing the unapproximated magnetic moment 
operator (|^) this overestimate is reduced to a factor ~ 2. Agreement with 
the empirical value is achieved only by taking into account the exchange cur¬ 
rent contribution. In the case of the transition 'ip{2S) —>■ rjc'j the impulse 
approximation underestimates the empirical decay width by a large factor. 
For this transition the exchange current contribution is the dominant one. 
The net calculated result falls only slightly above the upper limit of the un¬ 
certainty range of the present empirical result. There is no instance, where 
the non-relativistic impulse approximation result is close to either the empir¬ 
ical value, or to the net calculated result. This feature emphasizes the point 
that the Ml transitions are peculiarly sensitive to the strength and form of 
the conhning interaction through the associated exchange current. The ulti¬ 
mate reason for this sensitivity is of course the fact that there is destructive 
interference between the different contibutions to the transition amplitude. 
Those calculated widths, which are exceptionally small therefore also have 
the largest theoretical uncertainty, as even tiny modihcations in the model 
parameters could cause large relative changes. 


4.2 Spin-flip transitions in the bb system 

For radiative spin-flip transitions in the bb system the matrix element At in 
([T^ ) takes the form : 

M = -^{h + h}. (23) 

6mb 

The factor —l/6mb arises from the charge factor —1/3 and the factor l/2mb 
in the magnetic moment operator. The dimensionless integrals Ji and A are 
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defined as in eqs. ©• 


The calculated decay rates reveal that even in the case of the bb system 
the non-relativistic impulse approximation is unreliable, except in the case, 
where the initial and hnal states have the same degree of excitation. In those 
cases it leads to overestimates of the net calculated decay rates by factors 
~ 2. The exchange current contribution is large for all these spin-flip tran¬ 
sitions, and in fact the dominant contribution for transitions between states 
with different excitation number. 


Transition 

NRIA 

RIA 

RIA -|- conf 


23.4 eV 

21.2 eV 

18.9 eV 

23^1 ^ 2^So 

3.45 eV 
4.09 eV 

2.02 eV 
3.70 eV 

0.215 eV 
2.77 eV 

2^50 ^ 

0.74 eV 

2.25 eV 

0.0069 eV 

3^5i ^ BSo 
2^So 
3^Si 3^So 

2.82 eV 
1.41 eV 

2.06 eV 

3.56 eV 
0.35 eV 
1.83 eV 

0.126 eV 
1.09 eV 
1.32 eV 

3^So 

3^ So 2^Si 

1.15 eV 
0.296 eV 

3.50 eV 
0.529 eV 

0.347 eV 
0.253 eV 


Table 6: Decay widths for the bb system. The columns NRIA and RIA con¬ 
tain the results of the non-relativistic and relativistic impulse approximations 
respectively. The net calculated decay width obtained by combining the rel¬ 
ativistic impulse approximation and exchange current contributions is given 
in the column RIA-|-conf. 


4.3 Spin-flip transitions in the system 

The magnetic moment of the ^Si ground state of the system that is 
obtained from the expression is 


^ roo nl 

fj.{RIA) = -m„fXN / dpp^ / dz 
3 JO J-i 

/*oo rOQ 

drBp{r') 


dr r 
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\mh\Jl + vf 

3o{.'r'p)jQ{rp)ip{r). 



(24) 


Here (^(r) is the radial wave function for the ^Si state. If the velocity de¬ 
pendent “relativistic” correction terms are dropped from this expression it 
reduces to the static quark model result 




, rnp\ 

Q I “I fhV- 

6 I rric mb) 


(25) 


With the quark mass values in Table the static approximation value is 0.49 
Pat. When calculated from the unapproximated expression (p^ this result is 
reduced to 0.420 /xat. 

The exchange current contributions to this magnetic moment value is 
obtained from the spin-symmetric terms in the exchange magnetic moment 


expressions and for the conhning and gluon exchange interactions re¬ 
spectively. The explicit expressions for these corrections are 


.(C) = -^c 


p{G) = 


O ^ 1 9 

27 mcTrih //t^. 

Numerical evaluation of these two expressions give the results p{C) = —O.IO/iat 
and p{G) = O.OSO/xtv respectively. When these values are added to the im¬ 
pulse approximation result 0.420 pjq the net calculated magnetic moment 
becomes 0.35 /iat. The static quark model value therefore is expected to rep¬ 
resent an overestimate of the net magnetic moment value of about 23 %. This 
conclusion is in line with estimates of the exchange current and relativistic 
corrections to the static quark model predictions of the magnetic moments 
of the baryons, but in that case the situation is more complex because of 
the possibly substantial exchange current contributions associated with fla- 


mt 


+ 


+ 


mt 


+ 


mt 


drr^ip‘^{r)fiisf, (26) 


drrp‘^{r)pN- (27) 


vor dependent meson exchange interactions 15, 18 . 


The spin-flip transitions in the systems differ from those of the cc and 
bb systems in that they also obtain contributions from the gluon exchange 
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magnetic moment operator (0)- This contribution turns out in most cases 
to be weak in comparison with that from the exchange magnetic moment 
operator that is due to the conhning interaction This is shown in Table 
^ where the calculated decay widths for the spin-flip transitions between the 
(predicted) S'—states of the system are given with and without inclusion 
of the gluon exchange current contribution. The role of the gluon exchange 
current contribution is in every instance to increase the net exchange current 
contribution. 


In the non-relativistic impulse approximation we obtain spin-flip transi¬ 
tion decay widths that are close to those obtained in ref. [IT]. As shown by 
the comparison between the calculated decay widths that are obtained in 
the non-relativistic impulse approximation with the values obtained with the 
complete Dirac magnetic moment operator, the former approximation leads 
to overestimates by factors, which are typically larger than ~ 2. 


The role of the exchange current contribution is predicted to be large for 
the spin-flip transitions in the B'^ system. The net calculated result typically 
differs from that obtained in the relativistic impulse approximation by an or¬ 
der of magnitude, the only exceptions being the transitions between states 
with the same excitation number. 


Transition 

NRIA 

RIA 

RIA -|- conf 

RIA -l-conf -|- gluon 


138 eV 

91.4 eV 

42.5 eV 

51.4 eV 

2^^i ^ l^So 
2^Si 2^So 

62.1 eV 
43.0 eV 

375 eV 
27.0 eV 

5.40 eV 
3.08 eV 

46.1 eV 
3.94 eV 

2^So 

12.6 eV 

249 eV 

0.688 eV 

15.6 eV 

^ 21^0 

38.5 eV 

21.3 eV 

27.4 eV 

145 eV 
110 eV 
15.8 eV 

30.5 eV 
56.9 eV 
0.0145 eV 

129 eV 
85.7 eV 
0.0758 eV 

31^0 ^ 

31^0 ^ 

14.3 eV 
3.00 eV 

107 eV 
66.5 eV 

17.8 eV 

23.8 eV 

80.8 eV 

35.8 eV 


Table 7: Decay widths for the B'^ system 
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5 Discussion 


The present results suggest that the exchange current operator that is associ¬ 
ated with the scalar conhning interaction through current conservation gives 
a crucial contribution to the decay rates for radiative spin-flip transitions in 
heavy quarkonia. In most cases it interferes destructively with the matrix 
element of the single quark magnetic moment operator. As a consequence 
the value for these decay rates obtained in the impulse approximation differs 
from the net calculated decay rate by large factors. 

The results support the observation made in refs. that the rel¬ 

ativistic corrections to the single quark magnetic moment operator are of 
signihcant magnitude, although with the present unapproximated treatment 
of the relativistic magnetic moment operator, these relativistic corrections 
are moderated somewhat. 


The most important exchange current operator is that associated with the 
scalar conhning interaction. The gluon exchange current only contributes a 
small additional correction to the radiative decays of the Bc+ system. Both 
these exchange current operators involve excitation of intermediate virtual 
qq pairs. The exchange current corrections therefore do involve the negative 
energy components of the quark spinors. That such negative energy compo¬ 
nents are numerically signihcant has also been pointed out in ref. 0 , although 
in a quite diherent approach to the problem, based on an instantaneous ap¬ 
proximation to the Bethe-Salpeter equation. The present Schrodinger equa¬ 
tion approach with exchange currents, pioneered in refs. would appear 
somewhat more transparent, and as shown above appears to yield decay 
rates for the spin-hip transitions of heavy quarkonia, that compare more 
favorably with extant empirical values. The present Schrodinger equation 
approach may be derived from the Bethe-Salpeter equation by means of a 
three-dimensional quasipotential reduction in the adiabatic limit [|K . 


The present results suggest that the exchange current operator that is 
associated with the conhning interaction should play an important role in 
radiative transitions of heavy vector mesons to heavy pseudoscalar mesons, 
in particular the transitions D* D'f, D* Dg'y and B* B'y. 


Acknowledgments 


21 




DOR thanks Dr. Hakan Snellman for drawing attention to the Ml transi¬ 
tions in heavy qnarkonia. TL and CN thank the Finnish Society of Sciences 
and Letters for snpport. Research snpported in part by the Academy of Fin¬ 
land nnder contract 34081. 


References 

[1] S. Godfrey and N. Isgnr, Phys. Rev. D32 (1985) 189 

[2] E.J. Eichten and C. Qnigg, Phys. Rev. D52 (1995) 1726 

[3] G.S. Bah, K. Schilling and A. Wachter, eprint |iep-ph/9611226| 

[4] G.S. Bah, K. Schilling and A. Wachter, Phys. Rev. D56 (1997) 2566 

[5] G. T. H. Davies et ah, Phys. Lett. B345 (1996) 42 

[6] G.T.H. Davies et ah Phys. Rev. D56 (1997) 2755-2765 

[7] H. Grotch, D.A. Owen and K.J. Sebastian, Phys. Rev. D30 (1984) 1924 

[8] X. Zhang, K.J. Sebastian and H. Grotch, Phys. Rev. D44 (1991) 1606 


[9] 

[ 10 ] 

[ 11 ] 

[ 12 ] 

[13] 

[14] 

[15] 

[16] 


D30:1924, 1984 

J. Linde and H. Snellman, Nncl. Phys. A619 (1997) 346 |hep-ph/9703383 

D. Gromes Phys. Lett B202 (1988) 262 

F. Abe et ah (GDF Gollaboration) Phys. Rev. Lett. 81 (1998) 2432 

E. J. Eichten and G. Qnigg, Phys. Rev. D49 (1994) 5845 

W. Bnchmiiller and S.-H. H. Tye, Phys. Rev. D24 (1981) 132 

T. Barnes, Third Workshop on the Tan-Gharm Factory, J. Kirkby and 
R. Kirkby, Eds., Editions Frontieres, Gif-snr-Yvette (1994) 411 

K. Dannbom et ah Nncl. Phys. A616 (1996) 555 

J. L. Friar, Mesons in Nnclei, M. Rho and D. Wilkinson, Eds., North- 
Holland Pnbl. Go., Amsterdam (1979), Vol.2, 595 


22 





[17] C. Caso et al. (Particle Dat Group), Eur. Phys. J. C3 (1998) 1 

[18] C. Helminen, eprint Piep-ph/9711252| 

[19] F. Coester and D. O. Riska, Ann. Phys. 234 (1994) 511 


23 




